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We provide a detailed analysis of quantum field theory around a collapsing shell and discuss several 
conceptual issues related to the emission of radiation flux and formation of black holes. Explicit 
calculations are performed using a model for a collapsing shell which turns out to be analytically 
solvable. We use the insights gained in this model to draw reliable conclusions regarding more 
realistic models. We first show that any shell of mass M which collapses to a radius close to 
^— ^ , r — 2M will emit approximately thermal radiation for a period of time. In particular, a shell 

which collapses from some initial radius to a final radius 2M(1 — e 2 ) -1 (where e « 1) without 
forming a black hole, will emit thermal radiation during the period M < t < Mln(l/e 2 ). Later 
on (t ^> Mln(l/e 2 )), the flux from such a shell will decay to zero exponentially. We next study 
the effect of backreaction computed using the vacuum expectation value of the stress tensor on the 
collapse. We find that, in any realistic collapse scenario, the backreaction effects do not prevent the 
formation of the event horizon. The time at which the event horizon is formed is, of course, delayed 
due to the radiated flux — which decreases the mass of the shell — but this effect is not sufficient to 
' prevent horizon formation. We also clarify several conceptual issues and provide pedagogical details 

of the calculations in the Appendices to the paper. 
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(N| ■ I. INTRODUCTION, MOTIVATION AND SUMMARY 

> , 

Classical general relativity allows for solutions in which matter can collapse in a spherically symmetric manner, forming 
a black hole event horizon when viewed from the outside region. The collapsing matter hits a spacetime singularity 
in finite proper time thereby preventing us from using the classical equations to predict the future evolution of the 
system, as viewed by an observer collapsing with the matter. On the other hand, as far as any outside observer is 
concerned, the collapsing material takes infinite time (as measured by the stationary clocks outside) to reach the event 
' horizon and hence the formation of the singularity has no influence on the outside region. While this may appear 
, bizarre to the uninitiated, general relativists have learned to live with this dichotomous evolutionary behaviour of the 

■ system. 

. £h ! A new layer of complication arises when we go beyond the classical theory and study quantum fields evolving in 

■ the background geometry of collapsing matter. If a quantum field is in the vacuum state in the asymptotic past, then 
at late times an observer at spatial infinity will see a flux of energy from the collapsing matter corresponding to a 
thermal radiation with temperature Th = l/8irM jlj. (Throughout the paper, we will use Planckian units with G, h 
and c set to unity.) The thermal nature is closely related to the exponential redshift experienced by the modes when 
they travel from a region close to r = 2M to r — > oo and — to that extent — the thermal nature of the radiation 
depends on matter collapsing to a size r w 2M. This process has been extensively investigated and is well understood 
as long as we treat the quatum field as a test field (see Ref. for a pedagogical treatment and a review of the 
literature). 

But the radiation flux to r — > oo is 'real' in the sense that the observer can collect it and use it to do work (say, e.g., 
to heat up some water) . Energy conservation then requires that the mass of the collapsing body must decrease due to 
the outgoing radiation, since it is the only possible source of energy. In that case, we cannot treat the background as 
fixed but must solve for Einstein's equations with a semiclassical backreaction term added (which is usually modeled 
through the vacuum expectation value (VEV) of the stress tensor with the adjective 'vacuum' referring to the choice 
that, in the asymptotic past, the quantum field was in the vacuum state) and solve for the radiated flux simultaneously 
and self-consistently. Several authors have studied this phenomenon 0, 0, H, H, 0, B H Ell EH El> 0) EH EE El] and 



'Electronic address: asccm@tifr.res.in 
^Electronic address: paddy@iucaa.ernet.in 



2 



in particular we refer to the review by Brout et al. [3(. As long as the flux of particles, characterized by the VEV of 
the stress tensor is small and slowly varying (in Planckian units), the backreaction can be self-consistently accounted 
for and studied in the geometry of a black hole with a slowly varying mass. In such a case, the backreaction effects 
are not significant and can be incorporated in a perturbative matter. 

Recently there have been claims in the literature [17| that the number of particles reaching an observer at large 
distances from the black hole, at late times, can diverge. It has then been conjectured that the backreaction from 
this diverging number of particles might be sufficient to prevent the formation of the event horizon. These claims are 
clearly contradictory to the usually accepted view mentioned above, and other authors have also questioned this view 
(see, e.g. Refs. [H,[l9j]; see also Ref. Ho])- In the light of such claims, we wish to revisit this issue in this paper. 

Right at the outset, let us clarify one elementary — but potentially misleading — aspect of the problem. During 
the classical gravitational collapse, although the collapsing object will cross its Schwarzschild radius (and thus form 
an event horizon) in finite proper time, an observer at rest at large distances from the collapsing object will see the 
collapsing body approach its Schwarzschild radius only asymptotically as R(t) — 2M(1 + 0(e~*/ 4M )) where t is the 
usual Schwarzschild time coordinate of the outside metric 

dr 2 

ds 2 Schw = -(1 - 2M/r)dt 2 + (1 _ 2M/r) + r 2 dtf . (1) 

(Throughout the paper, we shall deal with situations in which a spherical body collapses without rotation.) In the 
absence of any backreaction, a calculation of the stress tensor VEV will show a steady flux at late times corresponding 
to a (nearly) thermal spectrum with a temperature Tjj oc l/M. If one insists on waiting for an arbitrarily long time 
without accounting for backreaction, we will clearly end up with an arbitrarily large number of particles reaching large 
distances at late enough times. Such a divergence, of course, has nothing to do with black holes, and will arise for 
any luminous object if we pretend that its mass does not reduce due to the energy it radiates. For example, if the 
sun could radiate a steady flux of photons for an infinite amount of time without changing its mass or constitution, 
it would emit an infinite amount of energy! That is, this situation will arise regardless of whether the event horizon 
forms or not and the divergence merely indicates the need to correctly account for the backreaction (or, more simply, 
energy conservation). 

In principle, it is therefore inconsistent to assume that the exterior geometry is described by an exact Schwarzschild 
metric for two reasons: First, the mass of the central object is changing and second, the presence of an outgoing flux 
even at large distances means that the geometry is not strictly asymptotically flat. One could now imagine a situation 
like the following : the mass loss of the collapsing object (in our case, the shell) implies that the radius to which the 
shell must collapse to form an event horizon, is shrinking. This means that the shell spends a larger amount of time 
outside the event horizon, all the while radiating particles to infinity. All this is likely to be further complicated by 
the fact that if the mass loss rate is high enough, the exterior geometry is likely to be very complicated. The radius 
which the shell is chasing keeps reducing, and the possibility arises that this runaway process ends with the shell 
completely evaporating before the event horizon is formed. If this situation is generic, it would mean that black holes 
typically do not form, contrary to the popular view. It is therefore worth investigating this scenario to settle the 
effect of backreaction, and in this paper we will do this using some simple models for collapse. The simplifications we 
introduce to make the problem tractable are the following: 

(a) We assume that the collapsing system is a thin spherical shell with internal stresses arranged in such a manner 
as to allow it to collapse along some given trajectory in the spacetime. Then the internal metric is Minkowski and the 
outside is Schwarszchild with matching conditions determining either metric in terms of the other and the trajectory 
of the shell p]]. 

(b) We will assume that all computations can be performed in the 1 + 1 sector of the spacetime which ignores the 
two angular coordinates. This corresponds to using only the s-wave mode of the scalar field and allows one to use the 
tools of conformal field theory to compute the VEV of the stress tensor. This assumption is unlikely to influence the 
conclusions of the paper and is justified in Ref. Q. 

(c) We will assume that the semiclassical backreaction can be modeled by the VEV of the stress tensor, < T a f, > 
with Einstein's equations modified to G a b = 8ir(T a b+ < T a b >). There is general agreement that such an equation 
should be valid in some suitable limit (though no one has rigorously proved it; but see Ref. [2l| ) and we will proceed 
hoping for the best. 

Given these assumptions, the problem can in principle be reduced to one of solving a set of equations. Given any 
trajectory of the shell, one can compute the VEV of T a {, everywhere using the conformal field theory technique. The 
problem then reduces to calculating (a) the backreaction on the shell trajectory and (b) the flux radiated to infinity. 
In practice, unfortunately, the equations turn out to be quite intractable and one needs to obtain insights into what 
is happening using simplified models. We do this along the following lines in this paper. 

We first consider the role of the event horizon in the emission of particles to large distances (ignoring the effect 
of the backreaction). We conclude that operationally, its role depends on when the observer at infinity sets out to 
detect the particles. To understand this result, consider two collapse scenarios (say, A and B) in which shells (each 
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of mass M) start from the same initial radius Rq with further evolution being different: (a) In Case A, the shell 
continuously collapses and forms an event horizon. (The shell crosses r = 2M at a finite proper time as shown by the 
clock on the surface of the shell but after infinite amount of coordinate time as measured by the clocks of the observer 
outside the shell.) Standard calculations indicate that at late times the asymptotic observer will see a thermal flux of 
radiation, (b) In Case B, the shell follows exactly the same trajectory as in Case A until it reaches close to r = 2M. 
Its trajectory then deviates from that of Case A; the shell progressively slows down, and asymptotically approaches a 
final radius r = 2M(1 — e 2 )™ 1 . Obviously no event horizon is formed in this Case B for any value of e, however small. 
It is also clear that at very late times the geometry is static in Case B with the shell staying at a fixed radius outside 
the event horizon and there should be no flux of radiation. But for an arbitrarily large period of time we can arrange 
matters such that the trajectory in Case B behaves similar to that in Case A when it is hovering just outside r = 2M 
(as seen by an outside observer). The outgoing modes which cross the shell and propagate to future null infinity will 
lead to a thermal flux in this case as well. In fact, this is an elementary consequence of causality. What an observer 
at infinity sees at some given event V (at time t b s ) can only depend on the behaviour of shell trajectory which is 
contained within the backward light cone of V ■ The future trajectory of the shell — in particular whether it settles 
down at r > 2M or goes on to collapse through r = 2M — should not affect observations at V '. 

We therefore expect the Case B to be characterized by three distinct phases. In the first phase, the shell collapses 
from a large radius to some radius close to 2M but larger than its final asymptotic radius. In this phase, the trajectory 
of the shell is identical to that of Case A. We, therefore, expect a flux of radiation being emitted to infinity which starts 
from zero and builds up to the thermal flux value characterized by the temperature Tjj oc (1/M) in the timescale 
of the collapse (which is 0(M) for nearly geodesic motion). In the second phase, the shell hovers close to r = 2M 
before settling down to its asymptotic radius. If we choose the trajectory such that the timescale governing this 
phase is O(M), then during this phase also the trajectory is similar to that of Case A, and now we would expect an 
approximately thermal flux being emitted by the shell just as though it is going to collapse to r = 2M. The crucial 
difference between the two cases arise in the third phase, during which the shell is asymptotically coming to rest. We 
will show by explicit calculation and a numerical analysis that during this phase the flux of radiation from the body 
dies down to zero exponentially. If the first two phases are characterized by a timescale 0(M), then the intermediate 
phase during which the shell emits approximately thermal radiation lasts for a timescale of the order of Mln(l/e 2 ). 

This explicit demonstration is one of the key results of this paper. This result shows that the thermal flux of 
radiation observed from a collapsing structure is not directly dependent on the formation of the event horizon. When 
the event horizon does form, the thermal radiation continues to escape for infinite amount of time if the backreaction 
effects are ignored. This is, of course, unrealistic since it would require an infinite source of energy. So in any realistic 
collapse scenario, the emission of thermal radiation is an idealization which ceases to be valid at sufficiently late 
times. Our model calculation shows that a similar effect can be mimicked by having a system collapse to a radius 
r = 2M(1 — e 2 ) -1 asymptotically. The approximately thermal radiation will emanate from the body during the time 
interval M <t < M ln(I/e 2 ). At very late times (t 3> Mln(l/e 2 )) the radiation dies down exponentially. 

The second key result of the paper is related to the computation of semiclassical backreaction and its effect on the 
formation of event horizon. We show that the effect of backreaction is essentially to delay the formation of event 
horizon by an amount Sv (where v is the in-going Eddington-Finkelstein null coordinate) given by 



where /?/, is the proper velocity with which the unperturbed trajectory crosses r = 2M and Lh is the luminosity 
of the thermal radiation. This result is intuitively clear : The delay in the event horizon formation is governed by 
the product of Lh and the timescale of the unperturbed collapse, which is essentially set by the initial mass of the 
collapsing object. Physically this is simply the naive bound one might place on the amount of mass that the collapsing 
object can radiate via semiclassical emission, during the collapse. The velocity at (unperturbed) horizon crossing /?/, is 
completely determined by the initial conditions (namely that the shell starts collapsing with zero velocity at a radius 
r = Rq at initial time). For example, for geodesic infall starting at r = Rq, one has /3 2 = 1 — 2M/Rq, and typically 
one expects [3h < 0(1) for more general trajectories as well. For all such trajectories the backreaction, which varies 
on a timescale ~ M 3 , cannot prevent the formation of the event horizon which occurs on a timescale ~ M; at best 
it can delay this formation by a time ~ 1/M. (One may wonder whether one can arrange matters for [3h to become 
sufficiently small to increase this timescale arbitrarily. This looks implausible and we provide detailed arguments in 
the relevant section eliminating this possibility.) 

It is possible to use these results to also put a bound on the total amount of mass that can be radiated away before 
the event horizon is formed and we find that it is given by M _1 lnM for most of the realistic trajectories. These 
results strongly suggest that the conventional wisdom as regards the formation of black hole is correct and that the 
semiclassical radiation does not prevent the formation of the event horizon. 

The plan of the paper is as follows. Section |TT] deals with the computation of the renormalized stress tensor VEV 
for a collapsing shell, in the absence of backreaction. We first recall some basic concepts of quantum fields in curved 
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spaces in Sec. Ill Al To set the stage, in Sec. Ill B I we analyse a trajectory which does not form a horizon and whose 
collapsing phase is null. While being unrealistic, this toy example admits of an exact calculation of the stress tensor 
VEV. In Sec. Ill CI we turn to the more realistic timelike trajectory discussed above, and analyse the behaviour of the 
stress tensor VEV in various phases of the collapse (which, again, does not form a horizon). 

In section IITII we include the effects of the backreaction of the stress tensor VEV, on the background geometry. We 
begin by recalling details of a self-consistent picture as presented by Brout et al. [3j, which describes the exterior 
geometry in the presence of a small and slowly varying emission from the collapsing object. We then present a simple 
calculation based on this picture, to derive the expression ([5]) for the delay in formation of the event horizon in the 
presence of backreaction. We argue that even for (classical) trajectories which slow down drastically as they approach 
the unperturbed Schwarzschild radius, it is highly implausible that backreaction can significantly delay the formation 
of the event horizon. Finally, we calculate an upper bound on the mass that can be radiated away due to semiclassical 
emission, using the extreme trajectories studied in Sec. [TXJ We briefly highlight our main results and conclude in 
section IIVI Appendix A contains pedagogical details of results from 2-dimensional conformal field theory which are 
relevant to our calculations, and Appendix B contains detailed proofs of various results that are quoted in the main 
text. 



II. FLUX OF RADIATION FROM A COLLAPSING SHELL 



We begin by recalling several results related to the quantization of a scalar field propagating on a fixed background 
geometry, which is taken to be that of a collapsing object. At the zeroth order we treat the scalar field as a test 
field so that the stress-energy tensor of the scalar field does not affect the background. For simplicity, throughout 
the paper we will consider the collapse of a thin shell of mass M, so that the exterior geometry is described by the 
Schwarzschild metric, while the interior is flat Minkowski spacetime. Since we will eventually be mainly interested 
in order of magnitude estimates and asymptotic behaviour of various quantities, we will focus on broad properties of 
the collapse trajectory without going into details of the shell stress-energy tensor, etc. 



A. Background geometry and mode analysis 

The exterior geometry of the shell is given by the Schwarzschild metric, which we write using the Eddington-Finkelstein 
coordinates, as 



2M N 
r 



d^cxt) = — (1 —)dudv + r 2 {u, v)dQ 2 



= -{l-2M/r)dv z + 2dvdr + r' l <Kl\ (3) 
where M is the (constant) mass of the shell and r is defined implicitly via 

e(V ~ U)/4M = ~ &r ' 2M 5 r>Rs > (4) 

where r = R s {t) is the trajectory of the shell with proper time r. The interior geometry is that of Minkowski 
spacetime, whose metric we write as 

ds( int) = -dT 2 + dr 2 + r 2 dfl 2 = -dUdV + r 2 {U, V)dVL 2 = -dV 2 + 2dVdr + r 2 dVL 2 , (5) 

where we have defined the light cone coordinates V = T + r, U = T — r, and we use the same r coordinate to label 
2-spheres in the interior and exterior. These metrics can be matched at the surface of the shell to get equations for 
U(u) and V(v) which are valid throughout the spacetime outside r = 2M. We will always define the shell's trajectory 
by the value of its r-coordinate as r — R Sl however we will treat R s as a function of cither u or v depending on 
convenience. The matching conditions are 

dV-dU = 2dR, ; dv-du= 2 '"'' : ' 



(1 - 2M/R S ) ' 

dUdV = (l - dudv , (6) 

with the differentials understood to be along the trajectory. For a timelike trajectory, these can be used to write the 
equation of the trajectory in the following alternative forms which will be useful later 

2R> s (l-U>) = U»-(l- 2 -f) , (7) 
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2R s (l-V)=(l-^j-V\ (8) 

where the prime and dot denote derivatives with respect to u and v respectively, along the trajectory. Once the 
trajectory is specified, these equations completely determine the metric both inside and outside and also allow us to 
relate the coordinates in the two regions. 

We next consider the quantization of the scalar field in this spherical geometry. Rigorously speaking, one needs to 
introduce the spherical harmonics Yi m (0, 0) and separate out the angular dependence of the modes. The resulting 
expressions will lead to a VEV of the stress tensor which cannot be evaluated analytically. To make progress, we 
shall concentrate only on the s-wave component of the scalar field and reduce the problem to one in the r — t plane. 
In such a 2-dimensional context, one can use the techniques of conformal field theory to evaluate the VEV of stress- 
tensor. We further ignore the effects of the Schwarzschild potential barrier which arises in the exterior even for the 
s-wave component. (See Refs. [H [H, [f| for a detailed discussion of the accuracy of these approximations, and Ref. 
for elementary arguments supporting the s-wave approximation.) The scalar field is therefore taken to satisfy the 
2-dimensional Klein-Gordon equation both inside and outside the shell, i.e., 

djjdvf = , (interior) , 

d u d v ip = , (exterior) . (9) 

We shall now briefly recall the procedure for quantizing the scalar field in the 1+1 spacetime and collect together the 
relevant formulas. The general solution of the Klein-Gordon equation in the exterior region is 

<p(u,v)=f(v)+£(u), (10) 

which after matching at the shell becomes 

<p(U, V) = f(v(V)) + e(u(C0) = f(V) + i(U) , (11) 

in the interior region. Imposing the reflection condition that the solution vanish at the center r = (V — U)/2 = 
then gives 

V = f(V) - f(U) = f(v(V)) - f(v(U)) , (12) 

in the interior, where v(U) is the function v(V) evaluated at V = U. Wc think of these solutions as wave-packets 
and study each frequency mode separately. Denote by (p™ modes which are positive frequency with respect to 
Schwarzschild time on past null infinity T~ . These modes will then define a vacuum (see below) which corresponds 
to the Minkowski vacuum on X~ . The relevant in- falling part of i/?™ is 

(onZ-), (13) 

1 VKUJ 

and is normalized with respect to the Klein-Gordon norm 

/oo 
dvft(id v )h, (14) 
-oo 

and similarly with u. The previous discussion on matching and reflection conditions shows that the outgoing it- 
dependent part of these modes, which reaches future null infinity X + , is given by 

Vl n 7= e "" G(U) > (onl+), (15) 

where the function G(u) is defined by the following chain : An ingoing mode labelled by v — v cn t behaves like 
~ e -^Ucnt m the exterior. It enters the shell at the event V en t, and now behaves like ~ e - luJv ont(v) m interior. The 
function v ent (V) is determined by the matching conditions at P en t- At reflection the mode becomes ~ e — *"«ent(f) ) 
and after exiting the shell at the event Pexit it becomes ~ e -»"»cnt(f(«)) i ^phe function G(u) is hence given by 

G(u) =v eat (U(u)). (16) 

The scalar field is quantized by defining a vacuum state | in) with respect to these modes by writing the field operator 

as 

/•OO 

<p= / ^(a^+h-c), (17) 
Jo 



G 



where "h.c." stands for Hermitian conjugate, and the are annihilation operators for the state |in), 

a u \in)= 0. (18) 

We will work throughout in the Heisenberg picture, so that the state does not evolve, and hence all expectation values 
must be computed in the state | in) . 

A key point is that on X + , the field operator has modes which are positive frequency with respect to G(u), not 
u. So the state |in) is the vacuum of the G- modes, or the "G- vacuum". Since G(u) is a nonlinear function of u in 
general, this G- vacuum will contain particles corresponding to the "u- vacuum" , defined as the state annihilated by 
modes which are positive frequency with respect to u. We will denote this vacuum by the state |out), annihilated by 
operators b\ corresponding to modes (p^ ut , such that on I + , the relevant outgoing part of <^ ut is 

^t~_L= e -^, (onI+), (19) 
V47rA 

and it should be clear that the corresponding in-falling part on T~ is 

^°"t 1 e -iAfc(iO on j- , (20) 

where if u = h(v) then v = G(u). The field operator ip can be expanded in the (/?^ ut modes as 

/>oo 

ip= / d\(bwT + h.c), (21) 

JO 

and a Bogolubov transformation relates the a w with the b\ . (For definitions and properties of Bogolubov transforma- 
tions see Ref. Q. We will not require these details.) To summarize, we have 

/>oo />oo 

ip= du (a^: + h.c. ) = / dX (6 A <^ ut + h.c. ) , (22) 
Jo Jo 

where 

in 1 / e-^" (onJ") 
Vu ~ X -e-^ G ^ (on 1+) 




(on J ), 

where, if v = G(u) then u = h(v), and we have G(u) = v cn t(U (u)) as discussed earlier. 

It is now possible to write down an expression for the time dependent part of the VEV of the stress tensor at any 
stage during the collapse. This is given by (see Appendix A for details of the notation and derivation) 

< r ->y<«> - s (f )'" d - (f y 1 ' 2 = ( r ->s--< r ->-. < 24 » 

which is implicitly trajectory dependent. It can be easily verified that in the limit when the surface of the collaps- 
ing object approaches its Schwarzschild radius, the function G(u) has the asymptotic form G{u) oc — 4Me~"/ 4A/ 
independent of the details of the trajectory. Hence 

du) du {du) ~64M 2 ' (2b) 

thereby leading to a flux which asymptotes to (7r/12)T^ at late stages of the collapse. For later use, we write 
( T uu )^ aj (u) in an alternative form. First note that since G(u) — v cn t(U (u)), differentiating using the chain rule gives 
G' = dG/du = (dv e nt/dV)(dU/du) = U'/V with the understanding that all functions of v such as V — dV/dv and its 
derivatives (see Eqn. (|2T|) below) are to be evaluated at v = v en t(U(u)). Using the identity 

F l,2 d 2 F -i,% = _1 [2^ 2 lnF- (^lnF) 2 ] , (26) 
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that is valid for any function F(x), where S 2 InF = d x [(d x F)/ F], and repeatedly using the chain rule, it is not hard 
to show that, 



' rp \ traj 
,-l-uu Iq 



(V) 



-1/2 



(27) 



v=v CIlt (U («)) , 



Due to the presence of the chain of functional dependences v cn t(U(u)), explicit calculations of the VEV will require 
knowing the functions v(V) and U (u). In principle, once a trajectory r = R s (u) or r = R s (v) is specified, the matching 
conditions Eqns. ([6]) can be solved to get U(u) and V(v). In practice, it is difficult to get both U(u) and V(v) in 
closed form, and the matter is further complicated by the inversion required to get v eni (V) at the point of entry V en t- 
The analysis can be done for a few special cases, and Brout et al. 0, [H| for example have an explicit calculation 
(in parametric form) of the time-dependent flux outside a collapsing star whose internal geometry is governed by a 
homogeneous dust (Ref. [3J, Appendix D). The asymptotic behaviour of attaining a constant flux is also displayed 
very nicely in that calculation (Ref. [3j|, Figure D.l). 



B. Trajectory with null in- falling phase without horizon formation 



Regardless of calculational difficulties, the general expression (|24|) highlights an important point : The presence of a 
flux of particles at large distances is governed entirely by the local dynamics (in outgoing time u) of the collapsing 
object. A flux will arise even in the situation where the object starts at some radius Ro and collapses to a final 
smaller radius Rf > 2M without forming a black hole, although the spectrum will not be thermal. Physically one 
would of course expect that the flux will only arise during the collapsing phase of the trajectory, and we will explore 
this feature next. We will first consider a situation in which a shell collapses from Rq to Rf along a null trajectory 
v = const., after which it remains static at r = Rf. Although this is a rather unrealistic situation, it does admit an 
exact calculation of the VEV of the stress tensor at any time. We will therefore present the calculations for this case 
first. Following this, in the next section, we will give results for a toy trajectory which remains timelike throughout 
and also does not form a horizon. In this case we will be able to provide order of magnitude estimates for the second 
term in Eqn. (|2T|) . while the first term will be exactly calculable. 

Recall that we are working in a geometry described by Eqns. (J3J)-(|H1) • Our example trajectory r = R s comprises 
three phases - 



u < 
< u < ui 

u > Ul 



R s (u) — Rq — const. 
v = const., V = const. 
R s (u) = 2M/(l-e 2 ) = 



const. 



with ui defined in terms of e as described below. We will ignore the derivative discontinuities at u = and u = Ui, 
since these can be smoothed away if needed, and in any case will not appear when we consider the more realistic 
timelike trajectory in Sec. Ill CI below. The parametrization of the final radius is chosen to ease comparison with this 
timelike case. For convenience we define the "tortoise" function x{r) as 



x(r) 



2Mln(r/2M- 1) ; r > 2M . 



(28) 



Using the trajectory equations ([6]), fixing some constants of integration and ensuring continuity at the transition 
events u — and u = u\, the functional form of the trajectory can be shown to be as follows: For u < (Phase (1)), 
we get 



2M 



1/2 



n // i = ( I ^1 j u - 2R ; V{v) = (l- ^f) 1 v - 2x(R ) (l-^f 



\ 1/2 

pi ■ - p i — i • i > ] ' 

Uq j \ n J \ jxq J 

For < u < ui (Phase (2)), we get 

v = 2x(Ro) ; V = ; U = -2R S ; u = -2x(R s ) + 2x{R ) , 
For u > ui (Phase (3)) we get 

[/(it) = e(u - ui) - 4M/(1 - e 2 ) ; V{v) = e(v - 2x(R )) ; iti = 2x(R ) - 2x{2M/{\ - e 2 )) . 



(29) 

(30) 
(31) 



The Penrose diagram for this trajectory is shown in Fig. [TJ There are five types of null rays (which we label I, II, ...V) 
relevant for the calculation of G{u) and hence the stress tensor VEV, characterized by the locations of entry and exit 
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FIG. 1: Penrose diagram for the trajectory described by Eqns. (|29|) - (|31[ ). We have set Ro = 12M and e 2 = 0.1. The thick solid 
line is the shell trajectory, and the dashed line is the timelike surface r = 2M. The interior of the shell is Minkowski spacetime 
while the exterior is Schwarzschild. The five types of null rays relevant for the stress tensor VEV calculation are also shown. 



events P e nt and V cx . These correspond to: (I) P en t and P ex in phase (1). (II) P cn t in phase (1), P x in phase (2). (Ill) 
Pent in phase (1), P ox in phase (3). (IV) The single ingoing ray v = 2x(Rq) or V = 0, which skims the trajectory and 
finally enters the shell at u — u±, with V cx in phase (3). (V) Pent and P ox in phase (3). For type (I) rays, we have 

2M\~ 1/2 , s / 2M N 1/2 



twC V 0= M__l V + const. ; U(u)=\l-—J u + const. , (32) 
and hence G(u) = u+ const, which leads to a vanishing flux. For type (II) rays, we have 

v cnt (V)=(l- 2 -^y 1/2 V + 2 X (R ) ; e -^M =e (H,-«o)/2M (^-2M) . Rs = _ u/2j (33) 



2M^ 1/2 



and hence 

G(u) = (l- U{u) + 2x{R Q ) , (34) 

with U(u) given implicitly by the last two equations in (|33|) . We will soon return to the flux obtained from this form 
of G(u). Note that type (II) rays all occur in < u < U\. For type (III) rays, we have 

/ 2Af\ -1 ^ 2 

Ventiy) = ( 1 ) V + const. ; U{u) = eu + const. , (35) 

\ R Q J 

which gives G(u) — e(l — 2A1/Rq)~ 1 ' 2 u+ const, which is also linear in u and hence gives a vanishing flux. Finally, 
type (IV) and (V) rays have 

v e nt(V) = e~ 1 V + const. ; U(u) = eu + const. , (36) 

also leading to G(u) = u+ const, and hence a vanishing flux. 

The only nonzero flux therefore arises for type (II) rays, in the interval < u < Ui. For (T uu ) tIlLS (u) defined in 
Eqn. (|2"4"1) , a straightforward calculation leads to 

< T ™> traj W = 7jrir^( 2 -^) ' °<"< u i> ( 37 ) 

487r R s {uy \ R s (u) J 
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Normalised flux 




FIG. 2: The behaviour of {T uu ) traj (w) normalized by the Hawking value (7r/12)T#, for the trajectory studied in Sec. Ill Bl We 
have set Ro = 12M. The solid and dotted lines show the flux for two values of e as labelled. The dashed line is the asymptotic 
Hawking flux corresponding to e = 0. 



and zero otherwise, where R s (u) in the given range is implicitly determined through 

£ -u/4M _ C (R„-R )/2M i R s ~ 2M ) _ (gg) 

(Rq — 2M) 

The behaviour of (T uu ) tra j(tt) normalized by the Hawking value (tt/12)T^-, is shown in Fig. ® for Rq = 12M and two 
values of e. The step-like rise and fall are due to the derivative discontinuities mentioned earlier, and do not occur in 
the timelike case which we will study next. The important feature to note is that for small but nonzero e, the flux 
attains the asymptotic Hawking value but eventually falls to zero. Eqns. |37|) and (|38|) show that the time taken to 
attain the Hawking value is governed by i?0i which is physically clear since the exponential redshift of the outgoing 
modes occurs only near r = 2M, so that for larger Rq, the shell spends more time radiating a smaller non-thermal 
flux. This argument can be easily verified numerically as well. 



C. Timelike trajectory without horizon formation 



We will now turn to a somewhat more realistic trajectory which is continuous, differentiable and timelike at all times, 
and which also does not form a horizon. The specific example we choose is a trajectory which remains fixed at r = Rq 
until u = 0, and asymptotically (as u — > oo) approaches a final radius Rf = 2M/(1 — e 2 ) > 2M. It turns out that this 
can be easily accomplished by appropriately choosing a function U' = dU jdu. [Since this calculation is for illustrative 
purposes, we will not worry about the shell dynamics needed to obtain the behaviour described below.] Consider then 
the following prescription for U' , 

U' = e+(^(l- 1 - e^j e- Q (") = e + Ae~ a ^ , (39) 

where < e < (1 — 2M/R ) 1 ^ 2 is a fixed constant, 

a(u) = 9 -^£duh(u) : (40) 

where 6(u) is the Heaviside step function and h(u) is chosen to have the following asymptotic behaviour 

h(u) = 0,u<0 ; h'(u) = 0,u< 0, 

h(u ->■ oo) -> 1 ; ft,'(tt->oo)->0, (41) 

and we require the asymptotic values for h and h! to be achieved exponentially fast, with a time scale determined 
by M. An example of a function h(u) which meets these requirements is h(u) — 9(u) tanh(K 2 u 2 ) for some constant 
k. The condition that the trajectory remain timelike throughout reduces to ~2R' S < U' or (1 — 2M/R S ) < U', which 
will in general impose a restriction on the allowed values of the timescale k for a given starting radius Rq- Physically, 
if the shell starts from rest at a larger radius Rq then it will take a longer time (at subluminal velocities throughout) 
to reach radii where the asymptotic exponential approach to the final radius begins. Hence a larger Rq will imply a 
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FIG. 3: Penrose diagram for the "asymptotic" trajectory described by Eqns. (|39fl - (|41[l . We have set Ro = 5M, e = 10 -2 and 
2Mk = 0.075. The thick solid line is the shell trajectory, and the dashed line is the timelike surface r = 2M. The interior 
of the shell is Minkowski spacetime while the exterior is Schwarzschild. The two null rays mark (by their entry points) the 
beginning of the infall phase and the approximate beginning of the asymptotic phase. 



smaller k. We have checked numerically that for h(u) = #(«) tanh(f« 2 u 2 ) with Rq significantly larger than 2M (say 
R > AM), the largest allowed value of K is K max ~ M~ 1 (M/R ) b with b fts 2 (and a very weak dependence on e), and 
in general we can expect b > 0. The Penrose diagram for such an "asymptotic" trajectory is shown in Fig. [3] 

The reason for this somewhat convoluted prescription is that we want to simplify the VEV calculations, which 
involve derivatives with respect to u, and it is therefore convenient to parametrize the trajectory using u. To see that 
the required behaviour for the trajectory R s (u) is reproduced by this U' , one analyses Eqn. which is reproduced 
below 



2^(1 - U') = U 
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1 



2M\ 



(42) 



The behaviour for u < is obtained correctly. For u > 0, consider the asymptotic regime (large u) where h ss 1 and 
h! w 0. In this regime U' 2 approaches e 2 from above, forcing (1 — 2M/R S ) to also approach e 2 asymptotically. In 
Appendix IB II we show that the asymptotic behaviour for the trajectory is 



1 



2M 

"r7 



■/(«) 



e -u/AM ^ £ 2 e - 



X/kM 



(43) 



where f(u) is an exponentially decaying function, the exact form of which depends on the value of e. Note that e is 
fixed and we are not taking ae^O limit. We then see that the required behaviour is being reproduced. 

As before, we wish to evaluate ( T uu )^ aj (u) at any given time along this trajectory. The first term in Eqn. (f2"T|) 
involving derivatives of U', can be easily computed exactly using Eqn. (|39[) and gives 



12tt 



{u'Y^dHu 1 )- 1 ' 2 



48tt 1 + {e/A)e a ^ 



1 - 2{e/A)e a ^ 
1 + {e/A)e ai - u ) 



2a" 



(44) 



It is interesting to note that in the case e = 0, i.e. for a standard trajectory which approaches r = 2M as u — > oo, this 
term is positive definite and gives rise to the Hawking flux (7r/12)T^, which can be easily checked by setting e = in 
Eqn. g3|). The behaviour of (8M) 2 [(C/') 1/2 9 2 (^')^ 1/2 ] f or h(u) = 9{u) tanh( K 2 u 2 ) with R = 5M , is shown in Fig. H 
Panel (a) shows curves for a fixed value of k with 2Mk = 0.075 (which ensures a timelike trajectory). The two solid 
lines correspond to e = 10~ 2 and e = 10~ 8 , while the dashed line corresponds to the standard case with e = 0, with 
the asymptotic value corresponding to the Hawking flux. In panel (b), the solid lines correspond to a fixed value of 
e = 10~ 8 , with 2Mk = 0.075 and 2Mn — 0.015. The e = curves for these values of k are also shown as dashed lines. 
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U-dependent backreaction eps=0 U-dependent backreaction 




u/2M 



(a)Varying e (b)Varying k 



FIG. 4: Behaviour of the [/-dependent backreaction term (8M) 2 [(U') 1/2 dl(U')- 1/2 }, for h{u) = 6{u) tanh( K 2 u 2 ), with R = 5M. 
Panel (a) : 2Mk = 0.075 is fixed. The thin solid curve has e = 10 -2 , the thick solid curve has e = 10 -8 , while the dashed curve 
is the standard case with e = 0. Panel (b) : For the solid curves e = 10~ 8 is fixed, while the dashed curves have e = 0. The 
thick curves are for 2Mk = 0.075 and the thin curves for 2Mk = 0.015. See text for discussion. 



We see that the initial behaviour of the [/-dependent term is practically identical to the standard e = case. 
Eventually the term involving a' in Eqn. (f44"l) becomes negative, leading to negative values for the right hand 
side at least for the chosen h(u). While this sign change might not be a generic feature, the late time behaviour 
(u 3> Kf~ 1 ln(l/e 2 ), see Eqn. (|43p) is generically seen to be an exponential decay ~ {e x / KM /eM 2 )e~ U / 4M , which 
directly follows from Eqn. (j44"|) . In Appendix IB 21 we discuss some additional features of the behaviour of this U - 
dependent term, which allow us to place a bound on the mass loss due to backreaction for this type of trajectory. 

We now turn to the second term in Eqn. (|27p . which involves derivatives of V , and an evaluation at V = v cn t(U (u)) . 
Unfortunately, unlike Eqn. (|44f> for the first term, it is not possible here to give an analytic expression valid at 
arbitrary times. The problem arises mainly because of the function v cn t(V) which — in general — has a complicated 
form. Recall that this function is determined by the matching conditions on the trajectory at the point of entry 
"Pent, and the evaluation is at V — U{u) where U(u) is determined by the exit point V cx n- Despite the lack of exact 
expressions, it turns out that we can make order of magnitude statements by realizing that we are essentially dealing 
with two kinds of ingoing rays characterized by the locations of V en t and V CX it on the trajectory : 

1- 'Pent occurs for u ont > but before the asymptotic phase so that u ont < k _1 ln(l/e 2 ). 

2. Vent (and hence Vexit) occurs in the asymptotic phase. 

We ignore rays which enter at u cnt < 0, since it is easy to show that v cnt (V) oc V in this case, giving V = const, 



and hence the second term will be exactly zero. For the choice of parameters used in Fig. 4(a) (with e = 10 -8 ), and 
fixing one constant so that Vt ra j(u = 0) = 0, numerically one finds that the ray which enters at w cnt = will exit at 
w ox « 14M, and hence the second term in Eqn. (|2T|) will vanish for u < 14M. In general if Vt ia j(u = 0) = then 
Mex(wont = 0) satisfies 



+ A e~ a ^du = 2R . (45) 



For rays of type (1), if u cnt < the requirement that the trajectory be well-behaved in the proper frame of the 
shell, allows us to argue that the second term in Eqn. (|2"7| must be small. For k _1 < u ent < k _1 ln(l/e 2 ), the same 
arguments show this term to be at most 0(M~ 2 ) and changing on a timescale 0(M) (see Appendix IB 31 for details). 
For type (2) rays we need to work entirely in the asymptotic regime, in which the function v eTL t(V) turns out to be 
approximately linear, and as we show in Appendix IB 41 the contribution from the second term of Eqn. (|27[) also decays 
exponentially. The details depend on the value of e, but the typical behaviour is 
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1 



,2 yrl 
v=v eDt (U(u)) e 1U 



e l/ K M e -Ku/4M . £ = 0(1) . e -«/4M <£ 2 e -l/«M ) (4 g } 



where u = u ox , the exit time. When e is significantly smaller than unity, this term will dominate over the first 
term (which falls like ~ (e 1 / KM /eM 2 )e _ "' 4M ) [22|]. This analysis demonstrates a result which one should intuitively 
expect on physical grounds : When the object is initially static, there is no outgoing flux. When the object starts to 
collapse, there is a flux of at most 0(M~ 2 ) which changes on a timescale O(M), and when the object approaches its 
asymptotically static configuration, this flux exponentially decays to zero. 
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III. THE SEMICLASSICAL BACKREACTION 



The presence of an outgoing flux of particles during the gravitational collapse of a body, immediately tells us that the 
assumption that the mass of this body is a constant, cannot be correct. The formal way to see this is to consider the 
energy conservation equation u b T^. a = at large distances from the collapsing object, assuming the Schwarzschild 
exterior (which is asymptotically flat). Taking the 4- velocity u b to be that of an observer at rest at fixed (r, 8, (j>), and 
integrating over a three volume of radius r, we have 

^ =jd 2 STl = WT[ = = -{T uu )p ee -L H . (47) 

The first equality follows from integrating the conservation equation, with d 2 S = r 2 sin 9d9d(j)\ the second follows 
from spherical symmetry; the third follows from assuming that the angular components of T Q {, vanish, so that the 
4-dimensional values of T tt , T rt and T rr are rescaled versions of their 2-dimensional counterparts; the fourth equality 
follows from replacing ^ 2d ^T a b by its semiclassical value computed in Appendix A. The last definition is for later ease 
of notation. 

In principle, this could lead to a situation in which the existence of an outgoing flux prevents the formation of the 
event horizon. The radiative flux decreases the mass of the collapsing body and thus decreases the effective radius 
of the event horizon which the collapsing surface is chasing. If the radius of the event horizon decreases fast enough, 
the collapsing surface will never be able to catch up with it and as a result the event horizon will not form. Our aim 
now is to analyse this situation more carefully and show that this scenario can be ruled out. Backreaction does not 
prevent the formation of the event horizon. 

In reality, what comes to our rescue is the following fact: the mass loss rate is set by the semiclassical backreaction 
which is governed by the ratio 1/M 2 in Planckian units, which makes it a very small number (e.g., for a solar mass 
object we have M Q ~ 10 38 ). It then turns out that under the assumption that the outgoing flux (and hence mass 
loss rate) is small and slowly varying, a straightforward ansatz is enough to estimate the effect of the backreaction 
on the background geometry self-consistently. Brout et al. [|[ present a calculation (which we sketch in Appendix 
IB 5|) which shows that the following late time scenario is self-consistent (see also Ref. [llj], and Refs. @, [l(| for the 
original papers): 

(a) The time dependence of the outgoing flux Ljj(u) at late times is determined solely by the time- dependent mass 

M{u), and is given by Lh ~ 1/M 2 , with u being an "Eddington-Finkelstein-like" outgoing coordinate. Lh is 
assumed to be small compared to unity (in Planckian units) and slowly varying. 

(b) The exterior geometry at large distances, in terms of the outgoing coordinate u, is the outgoing Vaidya solution 

given by 

d S 2 ^ tMlgcr = -^l-^^jdu 2 ~2dudr + r 2 dn 2 ; = -L H . (48) 

(c) The exterior geometry in terms of an ingoing coordinate v, at any distance, is approximately (i.e. up to terms of 

order O(Lh)) given by 

ds 2 cxt « - ( 1 - 2m( "' r) N ) dv 2 + 2dvdr + rW , (49) 



where the mass function m(v, r) is slowly varying in the entire exterior, in that 

£-<KL.) , £-<KL.), (50) 

with the transformation between the u and v coordinate being such that at large distances one recovers the 
outgoing Vaidya metric (|48|) with m(v,r) = M(u). (See Appendix IB 51 for details of the transformation.) 

(d) Self-consistency is demonstrated by showing that the VEV of the stress tensor (T uu ) tra ->(u) computed in this 
geometry does indeed behave as ~ M(u)~ 2 (l + 0(L H ))- 

The Brout et al. calculation (as acknowledged by those authors) is only valid in the regime where the flux Lh is 
small and slowly varying. There will inevitably be a phase in the collapse when the mass loss rate becomes significant 
enough that a perturbative expansion in Lh is no longer valid. A full fledged calculation of the backreaction in this 
regime, to our knowledge, has not yet been performed. Our interest however, is only to ask whether the backreaction 
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can delay the formation of the event horizon to the extent that it does not form at all. We will not worry about any 
significant backreaction effects that occur after the horizon has formed. 

It is important to note that all questions about event horizon formation must be asked in a reference frame where 
this formation occurs in a finite time in the unperturbed collapse. It is not possible to theoretically settle this issue if 
one insists on working entirely in the coordinates used by static observers at large distances, even though these may 
be the most natural coordinates to use, simply because even in the classical scenario, event horizon formation takes 
an infinite amount of time in these coordinates. 

Consider then a situation in which the unperturbed collapse trajectory of the shell crosses the radius r = 2M (with 
constant M) in finite proper time with a finite subluminal velocity. The exterior geometry for the unperturbed collapse 
is given by the Schwarzschild metric ([3]), while the interior is Minkowski spacetime (JSJ) - [It is also interesting to consider 
a trajectory whose in-falling phase is light-like, an extension as it were of the trajectory we studied in Sec. Ill Bl We 
have analysed such a trajectory in Appcndix lB 61 ] On parametrizing the timelikc trajectory using the shell proper time 
t (with ds 2 |traj = — dr 2 ), we find that the quantity dvjdr remains finite at horizon formation. It is then convenient 
to reparametrizc the trajectory using the ingoing Eddington-Finkclstcin coordinate v, as (r = R(v),V = V(v)). 
Assuming that horizon formation in this unperturbed case occurs at some time v = Vq, we have 

R(v)-2M = -k{v-v ) + O{{v-v ) 2 ) , as v ^ v , (51) 

where the constant k can be related to the proper velocity (3h = — {dR/dr)\^—2Mi as & = 2f3 2 (which follows from 
using the metric ijSJ) on the trajectory). The trajectory equation (jHJ) becomes 

2i?(l - V) = (1 - 2M/R{v)) - V 2 , (52) 

(where the dot is a derivative with respect to v along the trajectory) which shows that V also remains finite at 
R = 2M. 

Let us now incorporate the effects of a backreaction caused by a small outgoing flux Lh <C 1, and ask by what 
amount is the event horizon formation delayed. Specifically, let the trajectory now be (r = R(v),V = V(v)) where 
the interior metric is still Minkowski spacetime in (r, V) coordinates ([5]). The exterior metric, following Brout et a!., 
is given by Eqn. (|49|) . and the trajectory equation (|SJ) is easily shown to be replaced by 

2R(1 - V) = (1 - 2m(v, R(v))/R(v)) - V 2 . (53) 

For the exterior metric (|^|) . the event horizon (if it forms) is the last outgoing null ray r = r e h(v) and satisfies the 
outgoing null geodesic equation 

drch = 1 L _ 2m{v,r ch (v)) \ 
dv 2 \ r e h(v) J 

Note that this event horizon is distinct from the apparent horizon r — r a h(w) which is defined as the locus of events 
at which dr/dv — along outgoing null geodesies, so that r a h(v) = 2m(v 1 r a \ 1 (v)). However as we show in Appendix 
IB 5[ assuming that the flux Lh is small, we get 

2m{v, r ch (v)) = r ch {v)(l + 0{L H )) . (55) 

so that this distinction is irrelevant. We now use the smallness of Lh to make the following assumptions, which will 
turn out to be self-consistent at the end of the calculation, (a) We assume that the functional form of the trajectory 
is affected only by terms of order O(Lh), so that 

R(v) = R(v) (1 + 0{L H )) ; V(v) = V(v) (1 + 0(L H )) , (56) 

which is reasonable since the trajectory of the event horizon (r = 2M in the unperturbed case) is also affected by 
the same amount (see Eqns. (|55|l): and (b) the event horizon is formed (i.e. R(v) — r^v) is satisfied) at some 

finite v ~ vq such that vq = Vq + Sv where R(vq) = 2M and Sv is to be determined. We now evaluate the perturbed 
trajectory equation (|53p at v = Vq, linearize around v = vq assuming that Sv is "small", and use the unperturbed 
equation ([52^1 to obtain after straightforward algebra, 



°{^t Lh ) = °(m Lh ] * (:,T1 



where @h is the proper velocity with which the unperturbed trajectory crosses r = 2M and Lh is the luminosity 
of the thermal radiation. This result is intuitively clear : The delay in the event horizon formation is governed by 
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the product of Lh and the timescale of the unperturbed collapse, which is essentially set by the initial mass of the 
collapsing object. 

One might still argue however, that there could in principle exist unperturbed trajectories which slow down suffi- 
ciently while still outside r = 2M, that [3h becomes significantly smaller than unity. In such a case, it might be possible 
to delay the event horizon formation to such an extent that the backreaction itself starts changing significantly and 
becomes large, leading to a runaway process which might exclude the formation of the event horizon. While this 
situation cannot be excluded, the following argument does render it implausible. 

• First note that there is no semiclassical radiation in static geometries. The flux of particles arises when the 
geometry (characterized e.g. by the shell radius) changes by a significant amount during the time that an 
ingoing mode reflects at the center and exits the object. For this to occur the trajectory must be governed by a 
timescale linear in M, since this is the order of the time spent by the modes inside the object. For example in 
the "asymptotic" trajectory of Sec. Ill CI the maximum flux is achieved for times u k^ 1 when the timescale 
of the function a(u) is set by M. 

• Our calculations in Sec. HH and also e.g. the calculation in Brout et al.'s Appendix D, show that in any stage of 
the trajectory governed by the timescale M, the flux is expected to be at most 0(1/M 2 ), with the largest value 
expected only as the object approaches r = 2M. In Sec. Ill CI we also saw that in the extreme case when the 
trajectory is asymptotically slowed down to a halt, the flux in fact exponentially decays. One expects therefore 
that if a trajectory is slowed down by some smaller amount (i.e. not exponentially), the flux will have to lie 
between its maximum of 0(1/M 2 ), and zero. In particular, slowing down a trajectory cannot increase the flux. 

• In order to completely evaporate the collapsing object before the event horizon forms, what we need along the 
unperturbed trajectory is to sustain the maximum possible flux for the largest possible time. Naively one would 
want a flux 0(1/M 2 ) (which is the maximum possible) sustained for a time 0(M 3 ). The previous arguments 
show that this is unlikely to happen, since there is a trade-off between having a significant flux and remaining 
outside r = 2M for long enough. 



A. Bound on the total radiated mass 



For trajectories such as those discussed in Sec. HU and in fact for any trajectories that have such a form up to any 
finite time, we can place concrete bounds on the amount of mass that can be radiated away. We do this by calculating 
or estimating the integral AM = J (T uu ) tra 'du. For the trajectory of Sec. Ill Bl (T mi ) traj is non-zero only in the 
range < u < m, in which its integral can be explicitly performed to give 



AM -1 1 



M 48tt AM 2 



ln(e- 2 ) + ln(l - 2M/R ) + (2M/R ) - (1 - e 2 ) + | {(2M/R ) 2 - (1 - e 2 ) 2 } 



(58) 



We expect Rq/M to be significantly larger than unity. Also, since we are in the semiclassical domain, we can only 
expect to make meaningful statements for situations where the final asymptotic radius is not closer to r — 2M than 
one Planck length. This Planck length bound translates to a lower bound on e given by e 2 > 1/(2M + 1), so that the 
maximum mass that can be radiated away (ignoring numerical factors) is 

\AM max \ ~M- 1 In(M). (59) 

For the asymptotic trajectory of Sec. Ill CI the situation is trickier since we do not have exact expressions for the full 
backreaction, and also because the presence of the timescale k _1 causes a subtle interplay of competing effects. We 
show in Appendix IB 21 however, that even in this case there is an upper bound on the possible mass loss, which is in 
fact identical (up to numerical factors) to the expression (|59"]) . 



IV. CONCLUSIONS 



Since we have provided a detailed summary of results right at the beginning of the paper, we shall be brief in this 
section and will just review and stress the key conclusions which we have obtained. 

It is possible to explicitly compute the VEV of the stress tensor around a collapsing shell for any arbitrary trajectory 
of the shell. In principle, this calculation should be done in (1+3) dimensions but the results based on s-waves can 
be obtained by using conformal field theory techniques in the dimensionally reduced (1+1) case. Such an analysis 
shows that the time dependent geometry around a collapsing shell will lead to emission of a flux of particles to infinity 
though the spectrum, in general, will not be thermal. This result has nothing to do with black holes or event horizons. 
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When the shell does collapse to a radius close to the Schwarzschild radius, the relevant functions acquire the well 
known universal form characteristic of an exponential redshift of the outgoing modes. This exponential redshift leads 
to a thermal spectrum. In other words while any collapsing shell will lead to a radiation flux, such a flux will not, in 
general, have any universal characteristic. But if the collapsing shell goes close to its Schwarzschild radius, the flux 
becomes approximately thermal. 

If the collapsing shell follows a trajectory which does not form an event horizon and — instead — asymptotes to a 
radius 2M(1 — e 2 ) -1 , then we still get (approximately) thermal radiation during the period M < t < Mln(l/e 2 ). Of 
course, in the case e = 0, the collapse leads to an event horizon and we get the standard result that thermal radiation 
will be emitted at all late times corresponding to t > M. On the other hand, for e 7^ 0, at t > Mln(l/e 2 ) the radiation 
decays to zero exponentially as to be expected around a final static configuration. 

The above description does not take into account the effect of backreaction on the formation of event horizon. If 
radiation of energy can decrease the mass (and hence the effective radius of the event horizon) at a sufficiently rapid 
pace, then the collapsing shell might never catch up with the event horizon and a black hole may never form. This, 
however, does not happen essentially because the amount of radiation emitted during the collapse is not enough to 
decrease the effective radius of the event horizon sufficiently fast. This result is nontrivial and requires the careful 
computations which we have performed in this paper. 
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APPENDIX A: VEV OF STRESS TENSOR AND RELATED ISSUES 



In this appendix we present a pedagogical review of the calculation of the VEV of the stress tensor, using the tools 
of 1+1 conformal field theory, and highlight certain conceptual issues. The reader is referred to Refs. [1, U H, UH for 
details of the original calculations. The notation used here was laid out in Sec. Ill At 



1. Choice of vacuum state 



In the Heisenberg picture, there is only one vacuum state that one can work with, which is the |in) state. Nevertheless, 
it is common to find references in the literature to vacua defined with respect to various combinations of ingoing 
and outgoing null coordinates (for a review of these vacuum states see Ref. [23lj ) . One particular vacuum known 
as the Unruh vacuum takes on a special significance in the collapse situation. This is the vacuum defined to be 
positive frequency with respect to the ingoing Eddington-Finkelstein coordinate v and the outgoing Kruskal coordinate 



U = —4Me U / 4M 11]. As the surface of the object approaches its Schwarzschild radius, one finds that under the 



generic assumption that the trajectory in the (v,r) coordinates remains well-behaved at the Schwarzschild radius, 
we get G{u) -> -4Me- u/4M = U as u — ► 00. This behaviour arises due to the presence of an ingoing ray labelled 
by v = vh say, which after reflection is the last ray that can exit the object before the horizon forms (i.e. before 
the trajectory crosses r = 2M). The presence of this "last escaper" ingoing ray means that in the large u limit for 
the outgoing rays, the corresponding ingoing rays bunch together around v < vjj, and hence the function v on t(V) 
can be linearized. Additionally, assuming that the trajectory crosses the Schwarzschild radius at finite v = vq with 
a finite non-zero velocity dR s / dv\ VQ , it is straightforward to show using Eqns. that as u — > 00 (i.e. R s — > 2M), 
U(u) ~ e -u/4M _|_ (9( e -«/ 2M ). This shows that at late times G(u) cx e ~ M / 4M \ and in this limit, the G-vacuum is 
identical to the Unruh vacuum. 

The previous arguments provide a dynamical "reason" for "choosing" the Unruh vacuum at late times : the reason 
is simply the Heisenberg picture statement that the state does not evolve. The rest is taken care of by the dynamics 
and boundary conditions. 



2. Stress tensor VEVs : Regularization 



The object of interest is (in|T a (,|in), where T a b is the 2- dimensional stress tensor in the u — v plane. Provided the 
angular components of the stress tensor vanish, the four dimensional stress tensor can be obtained simply by rescaling 
the 2-d one by (47rr 2 ) -1 . The classical value of T a i> for a 2-d massless scalar field is d a ipdb(fi — (l/2)g a f,d l tpdiip, and 



1G 



it is then easy to see by expanding the field operator in the modes <^ ut , that (in|T a f,|in) corresponds ultimately to 
expectation values of the number operators bt b\ in the G- vacuum. Operationally one works with an expansion in the 
<p™ modes since the state |in) is the vacuum of the G-modes. 

However, the stress tensor VEV needs to be regularized, since even in Minkowski spacetime with \J Airunpu = 
e -iuv _ e — tum^ this object diverges. So what we want is 

< T ab ) rcn = (in|T afc |in) - "Minkowski value" . (Al) 



a. The accelerating mirror 



Let us first consider the toy example of a scalar field quantized in flat spacetime (ds 2 — —dudv + ...) with the boundary 
condition provided by an accelerating mirror which follows the trajectory v = G(u) for some G(u). Note that since 
there is no "interior" region in this problem, the complication of finding a function like v cn t(V) does not exist. The 
analysis is therefore simpler, and will help to set up some formalism which can then be easily adapted to the case of 
a genuine collapse situation. In particular one finds that the renormalized VEV required can ultimately be obtained 
in a fairly straightforward manner, without going into details of contour integrals etc. While the interesting case 
occurs when G(u) = — {\/a)e~ au for constant a, in which case reflection from the mirror leads to a thermal spectrum 
of particles in the G- vacuum, we will not require this form of G in the calculation below. The flat spacetime mode 
functions satisfying the reflection condition on v = G(u) are, 

f u = - 7 L= (e-^ v - e - lujG{u) ) . (A2) 

The "Minkowski value" for the stress tensor component T vv is simply (in.|T OT |in), while for T uu it is (out|T uu |out). 
This follows from the fact that the Minkowski vacuum has modes which are positive frequency with respect to u and 
v. This gives us 



"■" " ' 7 1>;| — (G') 1/2 d 2 (G')~ 1/2 , (A3) 



' T Y cn — • IT 
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where G 1 = dG/du. When G' = e~ au , the second expression can be obtained by explicitly writing out the mode 
expansions and performing the integrals involved (see Ref. [3J, Sec. 2.5.2). In general, the following regularisation 
procedure proves to be very useful, since it carries over to the gravitational case as well (see Ref. [|| for alternative 
regularisation procedures). 

One notes that the u-part of the Green function for the field in the it-vacuum |out), is (— 1/4tt) In \u — u'\ where 
one is considering two spacetime events (u, v) and (u',v'). Similarly, the Green function in the G-vacuum |in) is 
(— l/47r) In \G(u) — G(u')\. The stress tensor VEV in the state |in) can then be computed using a point split double 
derivative, as 

lim (±n\d u ip(u)d u > ip(u')\±n) 

u — >u' 

lim d u d u > (in\tp(u)ip(u')\in) 



= lim d u d u , [-J-ln|G(tt)-G(u , )| > ) , 

u~*u' \ 47T / 



(A4) 



and similarly (out|T uu |out), which gives the result in Eqn. (|A3|) after expanding the expressions to third order in 
(ti — vl). In order to get the result in a form which is generalizable to curved spacetime as well, one notes that the 
above argument for calculating the uu components can be inverted to get the GG components. One would do this by 
taking a point split double derivative with respect to G instead of u, and get 

(T GG Y cn = (in|T GG |in) - (out|T GG |out> - ~ Jim , d G G > (In |G - G'\ - In \u(G) - u(G')\) 



It is important to note that the u coordinate here corresponds to the Minkowski u. 
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b. Genuine gravitational fields 

Let us now turn to the original problem of a body collapsing to form a black hole. Due to mathematical similarities 
with the moving mirror problem, one can use similar ideas here as well, with some subtle differences 3, 8]. We still 
want a regularized version of (in|T a f,|in), the 2-dimensional stress tensor VEV. But now, there is no globally defined 
"Minkowski value" which one can subtract. Instead, one subtracts a locally defined object corresponding to the stress 
tensor VEV in a locally inertial vacuum. Given any spacetime with 

ds 2 = -C(G, v)dGdv + ... , (A6) 

such that we want expectation values with respect to the vacuum of the modes e~ luv and e~' lujG , one can show that 
locally inertial coordinates are given, at some event x = (Go, vq), by 

G c(g',vq) , . r cjGpyi 

Go C{G^) dG ; V = j V0 C(G^) dv • 

u and v are affine parameters along the radial null geodesies v = vq and G — Gq respectively, which pass through x, 
and the Christoffels in the (u, v) coordinates vanish at x 0]. At this event x one wants to subtract from (in|T a b|in), 
the quantity (l(x)|T a &|l(x)) where |I(x)) is the vacuum corresponding to the "locally inertial modes" e~ luJU , e~ luJV at 
x. 

The mathematics is now identical to the moving mirror case, and one finds 

1 / dii\ 1 ^ 2 „ 9 / du^ ' 



where the additional subscript on the VEVs reminds us that we are working in the G- vacuum (which is of course also 
the ^-vacuum). Eqn. (|A7j> can be used to remove all explicit reference to the locally inertial coordinates, and we find 
in general that 

(T G GYS n = -^C^d G C-^ ; { T vv Y£» = --±-CV*dtC-V\ (A9) 

1Z7T iZn 



For our case, the conformal factor is given in the exterior Schwarzschild geometry, by 

c-(i-?)* <w 

which follows from comparing the metrics in Eqns. ([3]) and (|A6|) . where u is the outgoing Eddington-Finkelstein 
coordinate. This simple prescription for constructing the stress tensor VEV is a consequence of considering only the 
s-wave contribution and further dropping the Schwarzschild potential barrier, which is what allowed us to approximate 
the scalar field modes in the exterior as propagating on a conformally flat geometry. See Ref. 0] for a discussion on 
how good this approximation is in practice. 



3. Changing vacua 

Let us see how the renormalized stress tensor VEV behaves under a change of the reference vacuum. This will allow 
us to write expressions which will clarify what happens in various physical situations. The structure of this object is, 
e.g. 



GG 



x rcn 
>G 



C 1/2 d%C- 1/2 - d G d G ' 



ln|G-G'|-ln|M 



(All) 



Choice of coords Choice of vacuum 



This shows that in any given vacuum, changing coordinates is trivial, one simply multiplies by the appropriate 
coordinate transformation factors. On the other hand, in a given set of coordinates, changing the vacuum is not as 
trivial, although it is straightforward. Let us suppose that we want to compute ( ) in the vacuum corresponding 
to some null coordinates (/, g) instead of (G,v), where we assume that / = /(G), g = g(v) so that the 2-d metric 
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remains conformally flat. [Later we will specialize to the case f = u, g = v.] Our previous results (Eqn. (|A9|) ~) tell us 
that 

{Tff) r f T = -^(C fg ) 1 / 2 d%C fg )- 1 / 2 ; (T gg }}™ = - 1 L(C fg )^d 2 g (C fg )-V\ (A12) 
with the subscript ( ) f g on the VEVs denoting the choice of vacuum, and where 

Cfg =C(G(f),v(g))^(f)^(g). (A13) 
In general for some function F(x), we have the identity 

F l / 2 dlF- 1 ' 2 = -I [ 2d 2 In F — (d x In F) 2 ] . (A14) 
Using this and the expression for Cf g in Eqn. (IA13|) . one can simplify Eqn. (|A12[) to obtain 

(^>™=(f) 2 (^>3--ii(f)'^?(fr i/2 

<^»^(|) 2 (^)--^(|)'^(|)"" 2 (A") 

In particular, for the "Boulware" vacuum defined through the modes e~' luJU and e -4 "", we have / = u, g = v. The 
Boulware vacuum attains physical relevance only in the static case when the star/shell does not evolve, and one has 
G{u) = u+ const (which can be shown using the trajectory and matching equations using R s = const). But the 
vacuum itself can of course always be defined even in the dynamical situation. For this vacuum, 

irp \ren _ / rp \ren \ ( dG\ ^ „ 2 f dG^ 



/ rp \ rcn j rp \ rcn 

\ 1 vv IB — \ 1 vv IG ■ 

As discussed earlier, in the limit when the object approaches r — 2M, we have G(u) 
hard to show that 



(A16) 

-4Me-"/ 4M , and it is not 



/ rp \ rcn 



' rp \ ren 
t J uti IU 



12 



-T 



H • 



(A17) 



where Th = 1/(8ttM) is the Hawking temperature of the collapsing object, and we have denoted the asymptotic 
G- vacuum by the subscript U (for Wnruh). Using the fact that the conformal factor for the Boulware vacuum is 
simply (1 — 2Af/r), where r is given implicitly in terms of u and v by the relation the stress tensor VEV in the 
Boulware vacuum can be explicitly computed as 



This gives us the relations 



,T UU )™(r) 



/ rp \ rcn 



1 



487T r 3 



^12 



3M 
r 



12 



T 2 



32M 3 48M 4 



en 
B 



' rp \ rcn 



(A18) 



(A19) 



= (T uu )T(r)-(T uu )T(r = 2M) 



12 



T 2 



12 



32M 3 



48Af 



41 



r I 



4M 
r 



12M 2 



(A20) 



where ( T u 



\ rcn 
IB 



\ rcn 
IB 



follows from the fact that d u r — —d v r. We see that (T, 



„,j 7 ^ n vanishes onr = 2M but 



attains the constant Hawking flux value as r — » 00. Also, since the coordinate transformation factor relevant to the 
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G —* U — > 0; r = const. < oo 
Static asymptotic 
observer at nnite r > 2M 


G^«^0;r^oo 
Static asymptotic observer 
at spatial infinity 


G -> U -> 0; (r - 2M) oc W 
Freely falling 
asymptotic observer 


/ rp \ ren 




n 
u 




/ rri \ ren 

\J-uu )b 


OO 


(if i . oo first) 


— oo 


j r-p \ ren 
{J-uu/U 


* 

Eon. (|A20|) 


* 

(tt/12)T£ 





/ rri \ ren 

{ J-uu )u 


OO 


oo 


* 

Finite, obs. dep. 



TABLE I: Various limiting cases of the stress tensor VEV in different vacua. 



sector is dU/du= (—U/4M) which goes to zero in the asymptotic limit, we see that the component (Iww)g n (r) — > oo 
in this asymptotic limit, for any finite value of r. [The limit hi — ► is distinct from the limits r — > 2M or r — * oo.] 

The object (Tuu)]j n is trickier. This object has a factor (AM /hi) 2 and a factor (1 — 2M/r) 2 in the asymptotic 
limit. Clearly it remains finite on the horizon if one takes the hi — > and r — ► 2Af limits simultaneously, assuming 
hi oc (r — 2M). Such a situation is relevant for an asymptotic freely falling observer maintaining approximately 
constant v, for example a geodesic observer who starts from rest at some very large distance. This result therefore 
shows that such an observer must see a finite flux at late times, although apparently this flux is not universal and 
will depend on the specific trajectory of the freely falling observer. 

For r = const. > 2M, the Kruskal coordinates are not the natural choice; such an observer uses the Eddington- 
Finkelstein (u,v) null coordinates. Hence, even though (Tuu)u n diverges as hi — > for fixed r, this is not physically 
relevant. The natural choice for such an observer in the hi — + limit is ( T uu )^ n (r), which as we see remains finite and 
attains the Hawking flux value at large distances r — > oo. Table [J summarizes these results. The boxes highlighted 
by asterisks show physically relevant limits. Not surprisingly, all of these are finite. 

4. Stress tensor VEV at arbitrary times 

The calculations in the previous section show that in the limit when the collapsing object approaches its Schwarzschild 
radius, the stress tensor VEV becomes time independent and takes the form given in Eqns. (|A20[) and (|A19p for the 
components ( T uu )g n and ( T vv )g n respectively, in the limit G — * hi. Physically this corresponds to a steady blackbody 
flux of particles at any radius r given by T[ = T vv — T uu = —(tt/12)Tjj. Showing that this flux is associated with 
a thermal spectrum requires a calculation of the Bogolubov coefficients relating the Boulware and Unruh vacua (see 
e.g. Ref. Q). We will not discuss this calculation since we wish to discuss backreaction, for which we need the flux 
and not the explicit form of the particle spectrum. In principle though, one already has the form of the stress tensor 
VEV at any stage during the collapse, since we have derived the relation (see Eqns. (|A16|) and (|A19|) ) 

<T uu )S n -<T„„)S» = — f^J d 2 J^j ={T m )%*(u), (A21) 

which defines the trajectory-dependent quantity (T uu )Q aJ (n) which asymptotes to (7r/12)Tj| at late stages of the 
collapse. Here (T vv )g n (r) = (T uu )g n (r) is given by Eqn. (|A"T8)) or Eqn. (|A"l9]l . 

APPENDIX B: PROOFS FOR VARIOUS RESULTS 

In this appendix we give calculational details for some of the results quoted in the main text. 

1. Asymptotic behaviour of the trajectory defined by Eqn. (|39[) 

As mentioned in the text, for U' given by Eqn. (|39|) , the trajectory equation (|42[) forces (1 — 2M/R S ) to approach e 2 
at late times. Let us assume therefore, that for large enough times the asymptotic behaviour is given by 

l-^ = e 2 + /(u) ; |/H|«6 2 . (Bl) 
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We will soon determine the condition to be satisfied by u in order for this to hold. Differentiating and rewriting Eqn. 
(ED leads to 



which are exact. The behaviour of U' in the asymptotic regime is well approximated by U' = e 2 + Ae^ u ^ M , where A 
is different from the A which appears in Eqn. (|39[) . The assumption that the asymptotic behaviour for the function 
h{u) in Eqn. (f^T|) is achieved exponentially on a timescale ac -1 , means that we can expect A = 0(e 1 ' K,M ). Recall 
from the main text that we have 

K< Kmax ~M-\M/R ) b ; b>0. (B3) 
Using the expressions in Eqn. (|B2|) in the trajectory equation (|4"2"| and linearising in / leads to 

4M/'(1 - e ) m (1 - e 2 ) 2 (2eAe~ u/iM - /) . (B4) 

Defining the variable s = u/AM , we have 

d s f + K 1 f = K 2 e- s , (B5) 

where 

(l-e 2 ) 2 

Kx = ^j, f- ; K 2 = 2eAX 1 . (B6) 

(l-e) 



f(u) = Ce~ Klu/4M + ^ 2 e -"/ 4M ; ifi^l, (B7) 
K 1 — 1 



Solving Eqn. (|B5|) gives 



with f = (C + K 2 u/ AM)e U / 4M when = 1, which we ignore since this solution forms a set of measure zero. The 
condition |/| <C e 2 shows that for consistency we must have 

u> K- x \n(l/e 2 ), (B8) 

which follows from the fact that for small e, Kij(K\ — 1) ~ A ~ e 1 /^ 1 . The constant C is not constrained by this 
analysis : it can only be fixed by evolving the full set of equations. The best we can say is that it must be at most 
of order ~ Ki/{K\ — 1), to be consistent with the condition |/| <C e 2 . Consequently the sign of / is also not fixed 
by this analysis. We see that f(u) exponentially decays and hence consistently reproduces (1 — 2M/R S ) — > e 2 at late 
times. 



2. Behaviour of the [/-dependent backreaction and bound on mass loss 

In this section we discuss some features of the backreaction term (f/') 1 / 2 9 2 ([7') _1 / 2 , and derive an upper bound on 
the expected mass loss due to backreaction in the "asymptotic" trajectory of Sec. Ill CI 



a. Analysing Fig. \J\ 

An important feature in Fig. [4] is the time (u = u* say) at which any given curve starts deviating from the e = case. 
This time depends not only on e but also on k. (For now we ignore the dependence on the starting radius Ro and 
assume k < K max .) We see that for a given k, smaller values of e will lead to a larger u*/M, which can be understood 
as the fact that it will take longer for a trajectory to "feel the effect" of a smaller e. 
On the other hand, for a given e, changing the value of k has three distinct effects. 

• Firstly, reducing k "stretches out" the standard e = behaviour of the backreaction, which is expected since 
the timescale over which the standard trajectory approaches r = 2M is governed by k^ 1 which is now larger. 

• Secondly, reducing k increases the value of u^/M at which the backreaction deviates from the standard case, 
which is also expected since a larger k^ 1 will slow down the rate at which a(u) increases and hence delay the 
time at which the effects of e 7^ are felt. 
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The third effect however is somewhat subtle and will play an important role below in determining a bound on 
the mass loss due to backreaction for such a trajectory. This effect is the fact that the increase in u^/M is 
generically less than the increase in n~ 1 /M. In other words, for two values k 2 < Ki 5 although we will have 
it*2 > we also see that k,2U*2 < Indeed, Fig. 4(a) shows that for the chosen values of K, K\u*i > 1 

whereas k 2 u^2 — 1. To understand the reason behind this behaviour, we note that the value of w* for some n 
and e is essentially determined by the condition e *-"*-* ~ A/e which follows from the structure of the right hand 
side in Eqn. (|44|) . Also, a(u) is an integral of a function h(u) which in turn is controlled by k. If we reduce k, 
then even though h(u) takes a longer time to rise to its maximum value of 1, the integral of h(u) can pick up 
a sizeable portion of its required budget of order ~ Mln(A/e) even for u < «2- This makes it plausible that 
although it*2 > u*i, we in fact obtain k 2 u^2 < It is difficult to give a more precise analytical argument 

for this observed fact. 



b. A bound on mass loss 



The upshot is that reducing n also reduces the maximum magnitude of the U -dependent backreaction term. We can 
use this fact to construct an upper bound on the mass loss due to backreaction along this "asymptotic" trajectory. 
Since the V-dependent term remains subdominant until u cnt ~ and thereafter is expected to follow a profile 

qualitatively similar to the [/-dependent term with slightly different timescales, for order of magnitude estimates 
we will not consider this term separately, and instead pretend that the [/-dependent term itself accounts for all the 
backreaction for the entire duration starting from u = 0. Our estimate for mass loss is a simple one : we ask for the 
maximum possible backreaction over the maximum possible time, and take the product. This will serve as a proxy 
for the integral f (T uu ) traj du (which is the mass loss at leading order, see Sec. Mil) . 

To begin with, we notice that for any value of k, the mass loss will be enhanced for smaller values of e, since 
reducing e will in principle allow the backreaction to be sustained at its maximum value of ~ M~ 2 (although the 



value of k may prevent this). This is clear from Fig. 4(a) So for any k, the largest mass loss will be achieved when e 
saturates its Planck scale bound (see main text, Sec. IIIIA"]) of e 2 > 1/2M. 

Having fixed e, we see that if we choose k appropriately so that the backreaction remains at its maximum value for 
times < u < ft ln(l/e 2 ), then the mass loss is approximated by AM ~ M ~ 2 kT x ln(l/e 2 ). The question now 
is what value of n will maximize this quantity. Naively one might think that making k smaller will help the case, 
since reducing k increases the time over which the backreaction is in effect. But as we saw in the previous subsection, 
this is not the whole story and there is a competing effect at play : reducing k will also mean reducing the peak 
magnitude of the backreaction, thus decreasing the mass loss. Clearly there is an optimum value for k which will 
maximize the mass loss. Fig. |4(b) and similar figures obtained by varying k for fixed e suggest that for values of 



1 max 



Rq significantly larger than 2M, this optimum value is in fact the maximum value K max , which is also borne out by 
numerically integrating the [/-dependent term for various values of k (and e as well). 

Finally, we note that K max is a decreasing function of the starting radius Rq, which we estimated earlier as n n 
M~ 1 (M/Ro) b for Rq significantly larger than 2M, where b > in general and 6«2 for h(u) = 0(u) tanh(« 2 u 2 ). If 
this scaling were true for all values of Ro, then a simple estimate for the largest possible value of the backreaction 
would come from setting k = K max with Rq/2M > 0(1) (for an e which saturates the Planck scale bound), i.e. by 
setting k ~ M up to numerical factors. However this scaling does not hold when Rq is close to 2M and K max in 
this case can actually be very large. Physically this is because it is possible for the trajectory to start its asymptotic 
phase very quickly without becoming superluminal. It might then seem that asymptotic trajectories which start close 
to r = 2M and very quickly enter their asymptotic phase will radiate away a negligible amount of mass. This is 
not true since the form of the backreaction in the initial stages of collapse also changes in this case : notice that the 
[/-dependent term (|44p depends on a" — h' /AM, and for large enough k, h approaches the Heaviside step function 
so that h! shows a spike close to u = of strength ~ M~ l . The bottomline is that even for small values of Rq for 
which K max can be very large, the mass loss is still of the same order of magnitude as it would be in our estimate 
M^kT 1 ln(l/e 2 ), with k ~ M _1 and e 2 ~ M _1 . The final upper bound on the mass loss (ignoring numerical factors) 
is given by 

AM max = 0{M~ X ln(M)) . (B9) 



3. Stress tensor VEV in initial phase of trajectory Eqn. (|39p 

Here we show that for type I rays which enter at times < u < k^ 1 ln(l/e 2 ), the stress tensor VEV is at most of 
order 0(1/ M 2 ) and changes on a timescale ~ M. The argument relies on simple order of magnitude estimates using 
the fact that the velocity in proper time j3 = dR s /dr in this phase satisfies j3 — 0(1) as we show below. [f3 actually 
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starts from zero, rises in magnitude, and then falls as the asymptotic phase begins, but for simplicity we will treat it 
as order unity throughout this phase.] We only need to concentrate on the second term in Eqn. (j2"T|) since we already 
have an exact expression (|44| for the first. 

Writing the trajectory matching equations §§§ in terms of the proper time dr 2 — — ds 2 | tra j and noting that /3 is 
negative for a collapsing trajectory, we have the following relations 

f = VTT]P + P; % = Jr ^ m (V^M/R^ + P), (BIO) 



^ = V^-^ i= {l -2M,R s) W^MjR^W-P). (BID 

Using these, one can write down expressions for U' and V in terms of (3 and R s . Further, in the initial phase of the 
trajectory we have U' = e + Ae- a{u ^ = 0(1), (1 - 2M/R S ) = 0(1) and hence dv/du = 0(1). This implies that (3 and 
hence V are of order unity. 

The second term in Eqn. (12"?]) can therefore be estimated as ~ Ty 2 where Ty is the timescale controlling V. This 
timescale is set entirely by the function U' due to the way we have chosen to define the trajectory. We can therefore 
estimate Ty 1 as 

Ty 1 ~ d u In U' ~ a ' = h{u) /AM , (B12) 

other factors being of order unity. Due to the required asymptotic behaviour of h(u), we can approximately take 
h{u) ~ for u < k^ 1 and h(u) ~ 1 for k^ 1 < u < k^ 1 ln(l/e 2 ). This is not too bad an approximation since, 
for example, if h(u) — tanh(«; 2 u 2 ), we will have h ~ k 2 u 2 for u <C k" 1 and h ~ 1 + 0(e~ 2K " ) for u > 
This estimate shows that the second term (J/'/V^) 2 ^ 1 / 2 ^ 2 ^ -1 / 2 is expected to be significant during the interval 
< M ont < ln(l/e 2 ) in which its value is at most 0(1/M 2 ) and changes on a timescale ~ M. 



4. Stress tensor VEV in asymptotic phase of trajectory Eqn. (|39[) 



Here we show that the second term in Eqn. (|27|) exponentially decays for the type II rays which enter and exit the 
shell in the asymptotic phase u 3> k" 1 ln(l/e 2 ) of the trajectory Eqn. (|39| . Using the late time solution (|B7j) . we can 
parametrically obtain v en t(V) by finding v(u)\-p Bllt and V(u)\-p ollt as follows : From Eqns. ^ we get 

d £ = 1 + T^mw, = 1 + ° (?^" ,4M ) - ■ (' + ° (k^' ,u ) ) • < B13 > 

where K, = 0(1) is determined from the solution (|B7[) and A = 0(e 1 ' K,M ) was defined in Appendix IB II above. This 
leads to 

V cnt (V) = ^ + (MAe-KV/AMe^ (B14) 

Further, since the exit point is also in the asymptotic phase, we have 

U = e(u + o(^e- u / 4M yj , (B15) 

which leads to 

v cnt (U(u)) = u + (Me^A . (B16) 

Using the expression for V from Eqn. (|B13|) above and replacing u (which is actually u Gn t) by v at the leading order, 
one can easily show that 

yl/2 d 2y-l/2 = Q f_A_ e -KW*M A f (B17) 
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and further, evaluation at v = v cn t(U(u)) simply replaces v by u (which is now u ex ). Also, at the leading order we 
have U' 2 /V 2 w 1, and hence the second term in Eqn. ([2T)) becomes 



^ v^dlv- 1 ' 2 
v 2 



V=Vent (U(u)) 



o 



A 



2 M 2 



-Ku/AM 



(B18) 



with the approximation valid for u 3> K 1 ln(l/e 2 ). 



5. The effect of backreaction on the background geometry 

In this section we sketch the arguments presented by Brout et al. which show that the late time backreaction is 
determined by the time dependent mass M{u) of the collapsing object. The reader is referred to Sec 3.4 of Ref. [H 
for further details. The argument proceeds as follows : The spherically symmetric exterior metric is written in the 
form 



ds^ t = -e 2V ^1 - 2m ^' r ^ dv 2 + 2e^dvdr + r 2 dn 2 , (B19) 
in terms of which the Einstein equations become 

where the right hand sides contain the two-dimensional stress tensor components. 

At large distances, the zeroth order calculation shows that only T uu — Lh survives. (We have already computed 
T vv which falls like ~ (M/r) 3 , and the energy-momentum conservation equation can be used to show that so does 
T uv .) Hence the metric must correspond to the outgoing Vaidya solution 

,9 ( 2M(u)\ , , , , 9 dM (u) 

^ext.lazger = ~ U ^ J ^ ~ ^udr + T 2 drf ; - -T uu , (B21) 

where we assume that Lh = 0(1/M 2 ) and varies slowly. Matching (|B19[) and (|B21|) at some large radius (say 
r > 0(6M)) gives the transformation between the (v,r) and (u,r) coordinates as 

2dr 

e^dv = du+ ; m(v,r)=M(u). (B22) 

1 — 2M(u)/r 

Using this to calculate the right hand sides of the Einstein equations (|B20|) in terms of Lh, one finds that at large r 
all three stress tensor components T r v , T v v and T rr , are of order 0(Lh) times some quantity of order unity. 

Near the apparent horizon r a h(v) = 2m(v, r a h(i>)), one uses the fact that r and v both behave like inertial light- 
like coordinates (see Eqn. (|B19[) ), to argue that T v v and T rr must be of order O(Lh) here as well. Integrating 
the conservation equation T r v r + T% v = from r a h to r > 0(6M) then gives T r v = O(Lh) near the horizon also. 
Specifically, one finds near the horizon (from Eqn. 3.86 of Ref. [|[) 

^ = -L H e+ + 0{ML H , v ) ; ^ = <D{L H ) ■ !tt = 0{L H /r), (B23) 
ov or or 

which shows that the metric function m varies slowly and also that tp can be safely set to zero in the calculation. 

The next step is to obtain the equation for outgoing null geodesies (ONGs). This is somewhat involved due to the 
presence of the apparent horizon where the ONGs change the sign of their "velocity" dr/dv. The analysis proceeds 
by first noting that as long as Lh is small and slowly varying, the location of the event horizon r = r c h(i>) is not very 
far removed from that of the apparent horizon. To see this, we write r e h(w) = r a h(^) + A(v). Differentiating with 
respect to v and using 2dr c \ l /dv = 1 — 2m(v, r e ^)/r e h (since the event horizon by definition is the last ONG to reach 
we find at leading order in Lh 

A = r ah r ahjt , - -1/M => r eh = r ah (l + 0(L H )) , (B24) 

since dr^/dv = 0(Lh)- A useful result which follows from this is 2m(i;, r e h) = r e h(l + O(Lh)), 

2m{v, r eh ) = 2m(v, r ah + A) = 2m(v, r ah ) + AO(L H ) + . . . = r ah (l + 0{L 2 H )) = r oh (l + 0{L H )) . (B25) 
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The ONGs are easier to analyse in terms of the coordinate x — r — r e b,(v), in terms of which the metric becomes 

(fe 2 xt = _ dv 2 2m (v,r ch + x)x {l + + + r2dn2 

r c h{r c h + x) 

The coefficient g vv vanishes on the event horizon (x = 0) in these coordinates rather than the apparent horizon, and 
the ONGs near the horizon resemble those in the unperturbed Schwarzschild geometry, namely we have for x <C r e h, 

/V 1 
—2m{v,r eh ){l + 0(L H )), (B27) 

along an ONG, and the following ansatz describes ONGs at arbitrary distances x, 

x f u du' 

5-2 — -2hxx + 5 = / —-+D. (B28) 

r c h(v) J m(vf) 

Here D is a constant of integration, and S — (D(Lh(Mx + x 2 )/M) which follows from integrating the ONG equation 
u = const using the form of the metric (|B26p . By requiring that at large distances the coordinate u be identical to 
the one appearing in the outgoing Vaidya metric (|B21[) . one finds rh(u) — M(u)(l +0{Lh)) using the transformation 
equation (|B22|) . 

This coordinate transformation is then used to calculate the quantity ( T uu ) lcn in the same way as was done in the 
case without backreaction, and to leading order one finds that the flux on the horizon vanishes quadratically and the 
flux reaching I + is (w/12)T^{u) where 

T H (u)~-^ = —{—(l + 0(L H )), (B29) 
8Trm(u) %ttM{u) 

which completes the argument (see Eqns. 3.91-3.95 of Ref. 0] for details). 



6. Delay in horizon formation for a null trajectory 

In this section we analyse the effect of backreaction on the horizon formation time for a null in-falling trajectory. Our 
classical trajectory in the absence of backreaction has two phases : 

1. u < : R s (u) = R — const. 

2. u > : v = const., V = const. 

The horizon is formed at finite (U, V). The trajectory solution is 

2M \ 1/2 



Phase (1): U{u) = ( 1 - — J u - 2R + AM ; 

2M\ 1/2 „ ,„ , / 2M X 1/2 



V(v)=^l-—j v - 2x{Ro) I 1 - —J +4M ; v = u + 2x(R ), (B30) 

Phase (2): V = AM ; v = 2x(R ) ; 

U=-2R S +4M ; u = -2x(R s ) + 2x{R ) , (B31) 

where x(r) = r + 2M ln(r/2M — 1) is the "tortoise" function and constants are chosen so that the horizon is formed at 
U = (with u — ► oo). The Penrose diagram for this spacetime is shown in Fig. [5] Causality implies that backreaction 
cannot have any effect on the spacetime region U < AM — 2i?o, i.e. while the trajectory is in phase (1). In phase 
(2) the trajectory will be modified and will no longer in general have V — AM — const. Say the modified trajectory 
(parametrized by U) is given by (V — V(U),r — R S (U)) in interior coordinates and (v = v(U), r = R S (U)) in exterior 
coordinates. The interior geometry is still Minkowski while the exterior is the Brout et al. approximation 

ds 2 xt rs - 1 ^— - dv 2 + 2dvdr + r 2 dtl 2 , (B32) 



where the mass function m(v,r) is slowly varying in the entire exterior. All this refers to phase (2) of the trajectory. 
The matching therefore leads to 

~ ~ I 2m (v R )\ 
dV -dU = 2dR s ; ~dU 2 - 2dUdR s = - ( 1 V sJ I dv 2 + 2dvdR s , (B33) 
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r=Rs r=6M 



FIG. 5: Penrose diagram (in the absence of backreaction) for a shell which stays at r = Ro = 3.5M until U = AM — 2Ro and 
then collapses along v — const., forming a horizon at U = 0. The shell trajectory and the r — 2M surface are labelled. The 
two null rays marking the onset of phase (2) and horizon formation, are also shown. 



and we will assume that the functional form of the trajectory is shifted by terms of order O(Ljj) so that 

R S (U) = R s (U)(l + 0(L H )) ■ R,(U) = --([/ - 4M ) , (B34) 

with M denoting the unperturbed mass of the shell. For the backreaction Lh, at the leading order we will take this 
to be given by the exact calculation of Sec. Ill Bl (see Eqn. (|3"7|l ). so that 

^-(T^W-jLf (2-2^). (B35) 

Although the calculation of Sec. Ill Bl was for a trajectory which does not form a horizon, the derivation can be easily 
modified to show that the stress tensor VEV (in the absence of backreaction) in the present case, will have the form 
(|B35|) at least for R s > 2M . Using this form for Lh is self-consistent so long as Lh <C 1, which is the case for all 
R s > 2M . We then have 

l + 2^=0(L H ), (B36) 



which gives us, at the leading order, 



% = 0(L B ) ; §=0(L H ). (B37) 



Now assume that the horizon forms at U = SU instead of U = 0, so that 

R S (5U) = r eh (v(5U)) ; v(SU) = v Q + Sv , (B38) 

where Vq = 2x(R a ) is the unperturbed value at which the horizon is formed. We would like to estimate 5U and Sv. 
Using Eqn. (|B37|) we have 

v(5U)=v + 0{L H )dU, (B39) 

J-2Ro+4Mo 

where the lower limit of integration is set by requiring that the trajectories in the presence and absence of backreaction 
be identical at the beginning of phase (2) which is at U = —2R + 4Mq. At this level of approximation, using the 



2G 



expression (IB35[) for Lh with R S (U) given by Eqn. (|B34[) . we find 

M (2M - U) 



SU »,r /r,».r rrN SU 



L H dU 



l-2R 0+ 4M 67T (4M - U f 

and hence, ignoring numerical factors, 



(B40) 

-2R +4M 



— = 0(M- 2 ) (1 + 0{SU/M) + O((M/R ) 2 )) . (B41) 

We can further assume r c h(t | ((5C/)) = 2Mq(1 + 0(M~ 2 )) since r e h(v) is a slowly varying function whose unperturbed 
form is the constant 2Mq, and since Lh = 0(M~ 2 ) around horizon formation. Then Eqn. (|B38[) gives us 

5U/M = 0{Af- 2 ) , (B42) 

so that \SU\/M <C 1 as expected. The signs of Sv and SU are not determined by this analysis. As in the timclikc 
case, the backreaction in this case cannot halt the horizon formation, but can only delay it by a time 0(1/M). 
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